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ABSTRACT 



With the standard electroweak interactions, the lowest-order coherent forward 
scattering amphtudes of neutrinos in a CP symmetric medium (such as the early 
universe) are zero, and the index of refraction of a propagating neutrino can only 
arise from the expansion of gauge boson propagators, from radiative corrections, 
and from new physics interactions. Motivated by nucleosynthesis constraints on 
a possible sterile neutrino (suggested by the solar neutrino deficit and a possible 
17 keV neutrino), we calculate the standard model contributions to the neutrino 
index of refraction in the early universe, focusing on the period when the tem- 
perature was of the order of a few MeV. We find sizable radiative corrections 
to the tree level result obtained by the expansion of the gauge boson propaga- 
tor. For z/g -|- e(e) ^ Ug + e(e) the leading log correction is about -|-10%, while 
for Ue + t'e(^'e) — * t'e + ^'e(^'e) the correction is about -|-20%. Depending on the 
family mixing (if any), effects from different family scattering can be dominated 
by radiative corrections. The result for z/ -|- 7 ^ z/ + 7 is zero at one-loop level, 
even if neutrinos are massive. The cancellation of infrared divergence in a coherent 
process is also discussed. 

PACS# 13.10. +q, 13.15.-f 



1. Introduction 

Neutrinos propagating in a medium could behave differently from in vacuum. 
The effect of the medium on the neutrino can conveniently be described by the 
index of refraction, n. Under certain conditions, matter effects could be so enor- 
mous that neutrino properties are modified drastically. The best known example 
is the Mikheyev-Smirnov-Wolfenstein (MSW) mechanism for the solar neutrino 
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deficit, ' in which the index of refraction of electron neutrinos propagating inside 
the sun could induce a sufficiently large effective mixing angle for them to oscillate 
into other types of neutrinos efficiently. 
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Neutrino oscillations could also be cosmologically important. One particularly 

5 
interesting example, pointed out by one of us (P.L.) some time ago, would be 

oscillations between a standard neutrino and a light sterile SU{2)i singlet neutrino, 

Vs- If certain conditions are satisfied prior to neutrino decoupling such oscillations 

can bring z/g into equilibrium and thus alter the expansion rate of the universe. 

In particular, for a wide range of masses and mixings one would have Ny = 4 

neutrino species in equilibrium, in violation of the bound Ni, < 3.3 (95% CL) 

from the observed abundance of ^He and other light elements. Subsequent more 
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detailed studies have reached similar conclusions. 

Theoretically, the possibility of having a standard neutrino oscillate to Ug has 
long been of interest. Many theoretical models based on the compactification 
of string theories and grand unification predict the existence of light sterile neu- 
trinos. A recent study shows that if the controversial 17-keV neutrino exists, 
laboratory, astrophysical, and cosmological bounds, unless significantly weakened, 
requires m,^^ = 17 keV and rrii,^ > rriu^. Thus, the MSW solution to the solar 
neutrino deficit would require Ve — ^s oscillations. It turns out that the mass and 



mixing parameters required for the MSW solution are close to the boundary that 
is excluded by primordial nucleosynthesis. Moreover, the exact location of that 
boundary depends sensitively on the index of refraction and other details. 

An accurate determination of the neutrino index of refraction is therefore an 
important step in the study of neutrino oscillations and their applications to solar 
neutrinos and cosmology. In the zero temperature limit, this can be done easily 
by a direct calculation of neutrino forward scattering amplitudes. The situation is 
more complicated if the thermal motion of scatterers becomes important. For that, 
lowest order calculations have appealed to a more complicated method derived 
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from finite temperature field theory. 

In this article, we wish to study the standard electroweak contributions to the 
neutrino index of refraction in the early universe. The epoch of interest is radiation 
dominated and has a temperature T of the order of a few MeV. At that time, 
according to the standard Big Bang cosmology, the universe consisted mainly of 
photons 7, electrons, positrons, and three types of neutrinos and anti-neutrinos. 

Motivated by the smallness of the baryon asymmetry, we will assume, for simplicity, 
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that in that era the universe was CP symmetric. As a consequence, there were 
equal amount of particles and antiparticles, with densities given by 



7- ^2 ^ ' 




4 




^ — N- - - - 


4^7, 



:i.ii 



where C(3) = E^=i "-"^ ~ 1-202. 

In terms of the standard electroweak interactions, the CP symmetry of the 
medium implies that the lowest order terms in neutrino coherent forward scatter- 



ing amplitudes cancel. This follows because the amplitude of, say, neutrino-electron 
scattering is equal in amplitude but opposite in sign to that of neutrino-positron 
scattering. As a result, the neutrino index of refraction is zero to lowest order. 
Indeed, previous calculations ( Ref. 11) have shown that the only nonzero result 
arises from the expansion of the W and Z boson propagators, and hence is sup- 
pressed by a very small factor of order q^/My^r, where q is the momentum carried 
by the W or Z. In the absence of family mixing, even such a tiny result vanishes 
for scatterings that take place between different families. 

The smallness of the tree level result calls for an investigation of its radiative 
corrections. In some places, particularly in resonance regions, theoretical solutions 
are sensitive to the value of neutrino index of refraction. Loop effects at zero 
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termperature have been discussed before. One finds that they introduce new 
interactions proportional to {a / n) (mj / M^^^) In^mj / M^^r) {i = e,fi,T) that break 
the universality of the neutral-current. In this paper, this effort will be carried over 
to the finite temperature case. Due to the aforementioned cancellation, radiative 
corrections now take a completely different form. 

The remainder of this paper is organized as follows. In the next section we re- 
formulate the lowest order calculation by introducing a method that uses standard 
field theory rather than the finite temperature field theory. Within this frame- 
work, higher order calculations are simplified to a standard calculation of radiative 
corrections to neutrino forward scattering amplitudes. In section 3, we present 
the leading log result for the z/g + e(e) -^ Ue + e(e) scattering. There, we begin 
with a discussion of the compatibility of coherent conditions with the cancellation 
of infrared divergence. Neutrino-neutrino scattering and its radiative corrections 
are discussed in section 4. Scattering between different generations of neutrinos 



and leptons is discussed in section 5. Neutrino-photon scattering is discussed in 
section 6. The imaginary part of neutrino index of refraction is given in section 
7. Our results are summarized in the final section and some technical details are 
presented in the appendix. The application of our results to cases of cosmological 
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and astrophysical interest will be discussed elsewhere. 



2. General Formula 



At zero temperature the index of refraction, n, is given by 



n-l = %Nf{0), (2.1) 

Po 



where po is the energy of the propagating neutrino, A^ is the number of scatterers 
per unit volume, and /(O) is the neutrino forward scattering amplitude. In the 
MSW formula, for example, /(O) arises from z/g — e scattering via the charged- 
current interaction, where in the rest frame of the electron /(O) = Gfpo/\^iv, 
and hence n = 1 + \/2NGf/pq. 

The fundamental formula Eq. (2.1) can be generalized to situations in which 
the medium has a finite temperature T. With the usual assumption that the 
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introduction of a finite temperature does not spoil the coherent condition, the 
only modification is to take the thermal average 



n-l = ^(Ar/(0)) = ^(iV)(/(0)), (2.2) 

Po Po 



where the last step follows from the usual assumption that scatterers are not cor- 



related. More explicitly 



3^ (2.3) 

(/(0)) = (iV)-^y"^n*(go,T)/(0), 

where 

is the occupation number density of the scatterer, /i is the chemical potential, and 
g* is the number of spin degrees of freedom. Within this simple scheme, physical 
quantities such as the forward scattering amplitude are calculated in the standard 
way (at zero temperature), and the final result is obtained by taking the thermal 
average according to Eq. (2.2) at the end. 
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This simple method makes the lowest order calculation of neutrino index 
of refraction trivial. As an illustration, consider the scattering V(,(j)) + e(p') -^ 
Ueip) + e(pO- The tree level Feynman graphs are shown in Fig. 1. Neglecting the 
electron and the neutrino mass, the scattering matrix element is given by 

/ -7(7 \ ^ — % 

Mo(z/ee -^ Uee) = - (^— =j _ ^2 _ j^fS \^^^V)laLUyXv)\\UvXv)l'^^Ue{jl)\ 

^9 \2 -i r_, , /x , ,. ^.2 



Ue{p')la{L - 2 sin ew)ue{p')][uySp)l°'LuyXp)]^ 



2cos%7 0-M|' 

(2.6) 

where dw is the weak mixing angle, and u^ and Uy^ are the electron and neutrino 
spinors, respectively. The first term in Eq. (2.6) arises from the -u-channel scat- 
tering (Fig. la), and the second comes from the t-channel (Fig. lb), where due to 
the coherent condition the momentum transfer must be zero. The -u-channel scat- 
tering matrix element has an additional minus sign, which can be understood as a 



consequence of fermion line crossing. One finds from Eq. (2.6) that the coherent 
forward scattering amphtude is 



m 



V2G 



fP-P 



27rp'o 



l + Asiii^Ow 2p.p' 



K 



+ ..., 



(2.7) 



where the first term in the parentheses is the same as that obtained from a local 
V — A interaction (i.e., {p — p') = 0), the second arises from the expansion of W 
propagator, and the ellipses refer to the higher order terms. It then follows from 

Eqs.(2.2) and (2.7) that 



(n-l) 



V2G 



F 



iy,,e^u,,e 



PO 



1 + 4 sin^ 9w 



A^. 



3M2, 



Po{Po)Ne 



:2.8) 



where Ng is the electron number density at temperature T and (pq) is the average 
electron energy. In reaching Eq. (2.8), we have assumed that the thermal motion 
of the electrons is isotropic so that {p') = 0. 

The matrix element for u^ip) + e{p') -^ J^eip) + e(pO (Fig. 2) can be obtained 
from Eq. (2.6) by changing Ug -^ Vg = Mg, p' -^ —p', and introducing an over all 
minus sign (due to the anticommutation of fermion fields). One finds 



(n-l) 



v^G 



F 



VeC^VeC 



PO 



1+4 sin^ Ow 



N, 



3M^ 



Po{Po)Ne 



(2.9). 



The sum of Eqs. (2.8) and (2.9) gives 



in-l)„^e{e)^i^,e{e) 



V2G 



F 



PO 



1 + 4 sin^ Ow 



{Ne - N,) - -^{{p'^)^Ne + {p'o)eNe) 



3M 



w 



(2.10) 



This is precisely the result obtained previously by the finite temperature field 
theory method (Ref. 11), in which one calculates one-loop contributions to neutrino 



self-energy in a background field. The present method is considerably simpler, since 
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it only involves a standard tree level graph calculation. 

Our approach also makes the underlying physics more intuitive. That the first 
term in Eq. (2.10) is given by the difference of the electron and positron densities 
is a consequence of the CP transformation property of vector currents 

UelaUe -* -VelaVe = -ullaU^ = -UelaUe- (2-11) 

Diagrammatically, from Fig.(l) and Fig. (2) the t-channel coherent scattering am- 
plitudes are equal in size but opposite in sign, and hence cancel completely. The 
same would also be true if the W^'s in Fig. (la) and Fig. (2a) were not carry- 
ing any momenta. In the neutrino-electron scattering the momentum carried by 
the W is p' — p, whereas in the neutrino-positron scattering it changes to p' + p. 
This relative sign difference compensates the over all sign difference discussed in 
Eq. (2.11), so that the second term of Eq. (2.10) becomes a sum rather than a 
difference. Finally, contributions from the longitudinal part of the gauge boson 
propagators are always negligible if the external fermion masses are small. These 
terms have the same structure as that in Eq. (2.10) times a multiplicative factor 
given by a ratio like rrig/M^. 

Similar observations can also be made for loop graphs that generate radiative 
corrections to the tree level result. Basically, all t-channel one-particle- reducible 
graphs will cancel, and hence we will omit them from now on. This includes all zero 
temperature correction results discussed in Ref. 14. Only the s- and -u-channel one- 
particle-reducible diagrams and the box diagrams may contribute to the radiative 
corrections. 



3. One-loop Corrections to the Ve + e(e) ^ z^g + e(e) 
Forward Scattering Amplitude 
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The order a/vr correction results are known. However, they have the same 

structure as the lowest order weak interaction, and thus the sum of the coherent 
amplitudes for neutrino-electron and neutrino-positron scattering cancels. 

Terms which may survive the cancellation must depend on [p ± p'Y where 
the two signs refer to neutrino-electron and neutrino-positron scattering, respec- 
tively. Comparing with the tree level result (~ Gp{p ± p'Y /M'y^) arising from 
the expansion of the gauge boson propagator, loop corrections are of the order 
{a/A'K)Gp{p±p')'^/M^r times a yet to be determined factor. Normally, this factor 
contains large log terms plus some constant terms of order unity. Since the ratio 
a/47r is already small, in the absence of logarithmic enhancement radiative correc- 
tions would be negligible. Thus, in what follows we will only keep those leading 
log terms. 

Explicit calculations (see below) show that there are only two scalar functions 
arising from loop integrals that have a large log 



1 

.2 /• ^2 I ^(^2 ^2\ I ^/^ i\^2 



^ / , ^ 9n y /" , / -, N , "^1 + ximi — mi) + x(x — l)q^ 
/i(mi,m2,M;g2) = A_y dxx{x-l)\n^ ^ ' ^^ ^ ^ 




(3.1) 
where mi^2 are two small masses, one for a neutrino and another for a charged 
lepton. In Eq. (3.1), M is either M\y or Mz, and q = p±p' is the momentum of 
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the virtual W or Z. The logarithmic enhancement arises because of the hierarchy: 



m 



2 I „2 1 ^^ ji #2 



l,2,lrl«M^ (3.2a) 



and hence ln(M'^ /mf 2)A^iM'^ /\p ± p'\'^) ^ 1. Consistent with our leading log 
approximation, terms proportional to I2 can be neglected, because the difference, 
which one must take in the end when summing over the neutrino-electron and the 
neutrino-positron scattering amplitudes, 

/2(mi, m2, M; (p + pf) - /2(mi, m2, M; {p - pf") 

f ml + x{m2 — ml) mf + x{m2 — ml) + x{x — l){p + p')"^ (3.26) 

J M^ m.f + x{ml—m.f)+x{x — l)(p — p')'^ 



does not have a large log. Effectively, this implies that we can neglect all small 
masses whenever possible in our calculation. A welcome consequence of such a 
simplification is that all graphs that contain scalars (in the Feynman gauge) can 
be neglected. As long as we keep only the leading log term, the result will be gauge 
invariant. We will assume that neutrino masses are completely negligible. 

3.1. Cancellation of Infrared Divergences 

A question which one might encounter in the study of radiative corrections is 

how do the infrared divergences cancel in a coherent process? In normal (inco- 

27 
herent) cases, the complete cancellation of infrared divergence is well understood. 

Although some individual Feynman graphs generate infrared divergent terms, they 

cancel when summed with the bremmstrahlung diagrams, which are also infrared 

divergent by themselves. 

11 



In a coherent process, however, bremsstrahlungs cannot be included. Thus, 
if a process is to remain as coherent at higher order, the cancellation of infrared 
divergences must take place among themselves. In the MSW formalism, the con- 
straint becomes even more restrictive. There, the neutral- and the charged-current 
interaction must be infrared finite seperately. 

The cancellation of infrared divergences for a coherent interaction has not 
yet received much attention, partly because there are not many places in which 
the physical process can be coherent. Marciano and Sirlin have investigated this 
question in the context of atomic parity violation. By explicit calculations, they 
show that the required cancellation does occur at one-loop level. In the following, 
we will show that this happens to all order in perturbation theory. 

Consider an electron scattering with an infrared finite but otherwise arbitrary 
potential and its QED corrections (Fig. 3). Suppose the forward scattering ampli- 
tude of the interaction is 

Mo{p,p) =Ue{p)To{p,p)ue{p), (3.3) 

where ro{p,p) is the interaction vertex, which is assumed to be renormalizable so 
that it is also free from ultra-violet divergences. 

At one-loop level, QED corrections introduce infrared divergences into some 
of the Feynman graphs. To isolate them, we recall that infrared divergences arise 
because photons are massless, and thus a slight acceleration of an external charged 
particle could result in an emission of an infinite number of soft photons. This 
implies that only those Feynman graphs in which a virtual photon line has both 
its ends attached to an external line (Fig. 3b to 3d) are infrared divergent (Ref. 
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27). The infrared singularity is generated by the photon pole 



P-V-T.^^-n6{e-\'), (3.4) 



where A is a ficticious photon mass. The one-loop QED corrections can now be 
written as 

6Mi{p,p) = 6M^ ip,p) + SM\^ {p,p) + 6M^ ip,p) + {infrared finite terms), 

(3.5) 
where (5M| {p,p) is the one-loop vertex correction (Fig. 3b) 

SMi'\p,p) = 

~^^^ I 77ru'^e{p)la —, TQ{p+k,p + k) — I'^Ueip) 



(27r)'* p.'-f + k.'-f — rrie ' p.'-f + k.'~f — rrig A;^ — A^ + ie' 

(3.6) 

and (5M| (p, p) and 5M^ {p,p) are the one-loop wave function renormalization 

contributions. 

In the on-shell subtraction scheme, one finds easily 

6M\^ {p,p) + 6M[ {p,p) = —Mo — In — ^ + {infrared finite terms). (3.7) 

^ /i 1 1 it:) 



"e 



The infrared divergence in the vertex correction can be separated out by considering 
the photon pole (Eq.(3.4)) contribution to the matrix element in the infrared region, 
k ^ 0, where by virtue of the finiteness of Mq{p,p) one finds lim.k_,QUe{p)To{p + 
k,p + k)ue{p) = Mq{p,p) and thus 

5MI '{p,p)\Q<ko<Loo = -e'^Mo{p,p) / 7^:34 ^^~r^^^(^^ " ^^) + {'infrared finite terms) 

2 

= —Mq(p,p) — In — ^ + (infrared finite terms), 
2n A^ 

(3.8) 

where ujq is an infrared cut-off often determined by the resolution of the experi- 
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mental apparatus. 

It is now evident that the infrared divergence in the wave function renormal- 
ization cancels the infrared divergence in the vertex correction. The net result is 
infrared finite 

2 

Mo{p,p)+ Yl SM^\p,p) = Moip,p)[l-^\n^)+..., (3.9) 

i—a,b,... ^ 

where the ellipses represent the other infrared finite terms. 

The factorization in the vertex correction (Eq. (3.8)) is crucial in reaching our 
final result. Now, by an elementary induction and by employing the same method 
presented above, one can readily show that at N-loop level, infrared divergences 
from the effective vertex (infrared finite at N-1 loop level) always cancel those from 
wave function renormalization. Thus, the forward scattering amplitude is infrared 
finite to all orders. This is to be expected, for if the electron forward scattering 
amplitude were not infrared finite, it would have implied that the electric charge 
is not a well defined quantity. 

3.2. QED Correction 

It is convenient to separate out the conventional photonic correction from the 

29 
weak interaction corrections. To do so, we employ the standard method to 

decompose the photon propagator 

-i Mir 



As explained in Ref. 29, the first term gives rise to a massive photon (7>) con- 
tribution. It's effect combined with those from W and Z will be refered to as the 
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weak interaction contributions. The second term is the same as for a massless 
photon (7<), but with an additional cutoff A = M\y as a regulator for ultra-violet 
divergences. The contribution of 7< (Fig. 4a) and (Fig. 4b) plus the photonic box 
(Fig. 4c) ( evaluated with the full photon propagator —i/k'^) gives the one-loop 
QED correction. The result for Ue + e ^ Ug + e is 



^mQ'^'' = ^£{[«e(p07A^^e(p')][«..(p)7'^^..(p)]^(3 In 

+ [Ue{p')l\l^Ue{p)][Uy{p)-i^Lu^{p)](l + 



(3.11) 
where q = p — p'. We have omitted all terms directly proportional to rrig. 

The result is finite in both the infrared and the ultra-violet region. Only the 
axial current (the second term in Eq. (3.11)) receives an order a correction. This 
is because QED does not introduce charge renormalization. It is known that the 
axial current term does not contribute to the neutrino coherent forward scattering 
if the electrons in the medium are not polarized. We will therefore neglect it from 
now on. 

A notable feature of (3.11) is that it contains a mass singular term InM^/mg. 

As we will see below, this mass singularity won't cancel even after we add the weak 

interaction corrections. This is not in conflict with the well known Kinoshita-Sirlin- 
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Lee-Naunberg (KSLN) theorem, however. The KSLN theorem follows because 
of the flniteness of Green's functions with non-exceptional Euclidean external mo- 
menta in a renormalizable theory. The non-exceptional momentum refers to con- 
flgurations in which no partial sum of momenta vanishes. Here, we clearly have 
an exceptional external momentum conflguration if one allows rrie -^ 0. The mass 
singular term is multiplied by g^ which vanishes in the limit m,^^ , rrie ^ in the 
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rest frame of the electron. Therefore, the KSLN theorem does not apply here. 
Comparing with the tree level result (taking only the first term of Eq. (2.6)) 



^0 = '-^J72^[^e{p'haLu,M]KMl''LUe{p')], (3.12) 



the leading log QED correction to the ratio q /M^y that determines the index of 
refraction (see Eq. (2.10)) is — (a/Gvr) ln(M^/'mg). The dominant contribution is 
provided by the wave function renormalization. Numerically, it is —0.9%. 



3.3. Weak Interaction Correction 



For reasons explained above, as far as the one-particle-reducible graphs are 

concerned only the -u-channel scattering graphs are of interest. This makes our 
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calculation very similar to that of /i-decay. Accordingly, we will employ the on- 
shell renormalization scheme of Sirlin (Ref. 31) 



cashew = j^. (3.13) 



For convenience, we orgainize our results into the self-energy, the vertex correction, 
and the box-diagram parts 



^M'^eak ^ ^^^^^^ ^ dMvertex + ^Mbox. (3.14) 

where 6Mseif and 6Mvertex also include the counterterm contributions and the 
wave function renormalization effect, respectively. 
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The self-energy and counterterm diagrams are shown in Fig. 5. The result for 
the sum of these graphs is known (Ref. 31): 



6Mself = Mo 



Aww{q^) - ReAwwiM^) 2(5e , c^ ^ fAzz{M^z) ^wwiM, 



f^i?e 



'W^ 



(3.15) 



where c = cosOw, s = sin^vK, and Mq is given by Eq. (3.12). The functions Aww 
and Azz are the coefficients of (^^jy in the WW, ZZ self-energies. These functions 

32 

have been calculated before. We will use the results given by Marciano and Sirlin. 

The vertex graphs are shown in Fig. 6. In evaluating the wave function renor- 
malization contributions we use the massive photon propagator 7> (the first term 
in Eq. (3.10)). The result is 



^ , ,'2 ,_ 1 M2 



6a+66+6c+6rf = -M^^{{- + T (1) - - - In ^j (^^-^^ 

- + r'(i)-i-in^)rl. 



(3.16) 
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To isolate the ultra-violet divergence, we have employed the dimensional-regularization I 
with ^l the ultra-violent cut-off and r(e/2) = 2/e + r'(l) + 0(e). 

Corrections to the vertex are obtained by evaluating the three-point functions 
from Fig. (6e) and Fig. (6f), in which one uses the full photon propagator. We 
find (details can be found in the appendix) 



«W^Mo^{(^r'(i,->„^)(i 



s2-c2 



2 4g2g2 

g2 r. M^ 



K 



)-^<'-iy^^'^-Ai 
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where 



/ 5 4c2/ 1 2c\ C^/ 8c2x 2 1 ,oic 



and 



xM^ + (1 - x)'^m'^^ 



/ 9n /" , / N, xMi + (1 ~ x) mi 

F{q^) = / dx2x{x - 1) In ^ ^ ^ ' ^ . 3.19) 

J x[x — l)q^ + (1 — x)m^ 



Combining Eqs. (3.16) and (3.17), we obtain 



«Ve«« =Mo^{Ul + r'(l) - i - to ^) + ^(5 + c') 1.-^ + 



s 



2 



-2 



2 ,.9 . (3-20) 

r 



^^ 



:^^^+(^-2^)^(^^)+^i} 



Comparing with the large log terms, one sees from Eq. (3.20) that the constant 
c' ~ 3 clearly can be ignored. Since we are interested in regions in which g^ jg 
space-like and Iq'^I ^ ml, in the leading log approximation one can make a further 
simplification 

Ftf)4ln-^. (3^21) 

Here, we have neglected the imaginary part that will be discussed in detail in 
section 7. 

There are a total of eight box diagrams (Fig. 7). A straightforward calculation 
shows that the leading log plus the order a terms are given by (details can be 
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found in the appendix) 



5Mbc 



ig a 



Ue{phaLUe{p)][Ury^{ph''LUu,{p)] 



X 



+ 



5c^-3s4, 2 15-24s2c2 ll + 20s2 + 56/ q'^ , q 

\nc^ -\ :-^^^ \ :r-^ ^^ In 



2s4 



is'^c'^ 



24s^ 



M^ Ml, 

l^^^^[Ue{p)laLUe{p')][UyXp)lpLUyXp)]<f(l^ ^"^ ^^ 



'W 



I2s' 



(3.22) 



In reaching this simple resuh, we have made the approximation 



h{mj,^,me,Mw]q 



1 q^ ^^ q^ 



6M2, M2,' 



(3.23) 



and ignored all non-leading log terms. The box diagrams also generate terms which 
do not contribute to the neutrino coherent forward scattering amplitudes. They are 
of the form [ue'^al^Ue][up^'^"'LupJ\ and q"'q^[ue'^al^Ue][ui,^'^pLupJ\ and have been 
omitted from Eq. (3.22). The result can be written in a more compact form by 
using the relation 



q'^q^[UelaUe][Uv,ll3LUj,J\ = —[Ue'JaUe][Uu,l'^ Lu„^]. 



(3.24) 



It then follows that 



SM 



Box 



'Ue{p'haLUe{p')][U:,^{p)-f'^LUrySp)] 



2M^47r 



5c2-3s\ 2 15-24s2c2 2 ,^ 2 . 4n 9 , 9 
' T^ + Tr^ 1 + 5 +4s^hVlnT79 



2s4 



In c^ + 



6S^C^ 



<-w 



+ ..., 



w 
(3.25) 



where the ellipses represent those (incoherent) terms of no interest to us. 
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Among the eight box diagrams, the first four also appear in /x-decay (with 
appropriate change of external lines). The order a correction of these graphs is 
given by the first term of Eq. (3.25). From Eqs. (3.20) and (3.25), one can see that 
the coefficients of g^/M^ in ^Mv^eriea; and in 5Mbox are almost equal but with the 
opposite sign. This reults in a large cancellation between these two contributions. 

3.4. Combination of Results 

The combination of all weak interaction corrections plus the tree level result 
(Eq. (3.12)) can be expressed in terms of the weak coupling constant G^ = 
(1.16639 ± 0.00002) x 10^^ GeV^'^ determined very accurately from /x-decay by 
the substitution (Ref.31) 



1 + Ar) = ^, (3.26) 



where 



^^ _ ReAww{M^) - AwwjO) _^^c^ (Azz{mI) Aww{M^, 






M2, e sH M| M^ 

a r4 /2 „,,., , Ml \ /7 „ 9\lnc2 






(3.27) 



It then follows from Eqs. (3.12), (3.14), (3.20), (3.25) and (3.26) that 



r 0-/15 — 'yAr^c^ 

Mo + (5M^^«^ = -^2v^G^[«e(p')7a^«e(p )][^..(P)7"^«..(P)]|1 - — ( 

+ 



47r V Ss^c^ 



—^ 1 + Ai?g^ +^( o 9 ln^-21n— ^ 

M^ L ^^ ^ 47r V 3s2 g2 ^2 

(3.28) 
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where 

AR( 2) - ^w ( ReAwwJM^) - Awwjq'^) _ ReAwwjM^) - Aww{0) \ 
^^ ^ - g2 Iv M2, - g2 M^ J 

- q2 + M^ +^(^^VJ 

(3.29) 

is the W self-energy contribution, with the understanding that only the real part 

of the result will be kept. 

The residual order a correction (the second term in Eq. (3.28)) arises because 
we have omitted all the t-channel one-particle-reducible graphs. Also, the last four 
box diagrams do not appear in the calculation of /x-decay. In any case, this term 
does not contribute to the index of refraction. 

Terms in the squared brackets will contribute to the neutrino index of refrac- 
tion. The first term is the tree level result obtained by the expansion of the W 
propagator. The second term, AR{q'^), is due to W self-energy and counterterms, 
and the rest is the combination of vertex and box diagram corrections. 

The large cancellation between the vertex and the box diagrams makes their 
contribution rather small. The sum is about +2% for |g2| = i MeV'^ and it 
becomes negligibly small, +0.1%, for |g2| = 400 MeV'^. 

The largest contribution comes from W self-energy. For convenience, we will 
decompose it into the bosonic, the hadronic and the leptonic parts 

AR{q^) = Ai?(^)(g2) + Ai?W(g2) + Ai?(^)(g2). (3.30) 

The bosonic contribution is negligible, because it does not have any large log term. 
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Explicitly, from Eq. (3.29) and the given self-energy function in Ref. 32 we find 



1 



^^•"('^)-s(i-5?^?/''^^(^-«0. p-31) 







where 



F{x,0 






2c2 2c2 2 



.2^2 



2s^ - 16c^ - lOc^ - 1 1 + 4c* + 16c^ 20c^ + 1 2^ ^ c^^;^ -x + 1 



In- 



c^a; — a; + 1 



2 — ,^ + ,^a; — a;2 x"^ + ^{1 — x) x{l — x) ^ s^x{l — x) 



X + ^{1 — x) X + C,{^ — x) c^x'^ + {1 — x)' 

(3.32) 

and ^ = m^i/My^. Here m^ is the Higgs mass. 

The hadronic correction seems to provide the biggest contribution, but its cal- 
culation is complicated by strong interaction physics. The conventional wisdom 
here is to employ dispersion relations to relate A^ \q ) to some physically mea- 
surable quantities. While a detailed study of this problem lies beyond the scope of 
the present paper, we will estimate the result by employing the constituent quark 
mass. We find for |g2| ^ rnf (here rrii is a constituent quark mass) 



1 



mfx + m^(l — x) — My^x{l — x) 

X{1 — X) In -K 7T- 

^ - miX + mAl — x) 

«J ^ ]^ ' 

mfx + m3{l—x) mfx + w?Al — x) — My^x{l — x) l 

2M2^ ^ mfx + m^-{l-x) 12. 

(3.33) 

where Vij is the KM matrix element. The sum is over all three families of quarks 

with iriu^d ~ 300 MeV , iris ~ 450 MeV , etc. The second term in Eq. (3.33) is 

obviously negligible for the first two generations. We present it here is to show 

that it is also negligible for a heavy quark. Indeed, for mt ^ Mw, it becomes 
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x{l — x)/2, which is much smaller than the first leading log term. Numerically, we 
find that ^R^^\q^) is about +5%. 

The calculation of the leptonic contribution is much less ambiguous. We find 
that the leading log contribution is (neglecting possible mixings in the lepton sec- 
tor) 

A«'"ta^) - ^ E i^e/...(l-.)l. "f^g^' . ,3.34) 

Numerically, it varies from +4% to +3% in the region 1 MeV"^ < |g^| < 400 MeV"^. 

Now, the sum of all weak corrections plus a —1% QED correction varies from 
about +10% to about +7% in the region of interest discussed above. Our results 
are summarized in Table 1. 

The result for z/g + e ^ z/g + e can be obtained from that of neutrino-electron 
scattering by changing Ue'JaLue to —Ue'JaRue and q -^ —q (valid only if the 
neutrino and the lepton mass is negligible). The outcome is that the constant 
term (g^ independent) has the opposite sign, whereas the g^/M-^ term remains 
the same. 

To obtain the index of refraction, we now need to take the thermal average 
of Eq. (3.28). Since the variation of radiative corrections in the region of g^ of 
interest is not significant, in the leading log approximation we can write the index 
of refraction in a A^g = -^e niedium as 

(n - l)uMe)^Ueeie) = -rT72"'^M^ePo(Po) (1 + ^u^e), (3.35) 

where 6i^^e is simply given by the non-averaged result discussed above; it is between 
about +10% to about +7% (see Table 2 for detail). Results for other values of g^ 
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can easily be obtained from the analytic formulae given above and in the appendix, 
provided the leading log approximation can still be justified. 



4. Neutrino- (anti) neutrino Scattering 

Fig. 8 shows the tree level Feynman graphs for a diagonal (in family space) 
neutrino-neutrino and a neutrino-anti-neutrino scattering. Fig. 8a and 8c are the 
standard t-channel forward scattering graphs induced by the neutral current inter- 
action, where the coherent condition requires that the virtual Z cannot carry any 
momentum. This feature is shared by all t-channel one-particle-reducible graphs. 
As explained before, the sum of coherent amplitudes of these t-channel graphs can- 
cels for V ^ V ^* V ^ V plus V ^ V -^ V ^ V. Therefore, we will ignore them from 
now on. 

Due to the identical particle property, there is another angle for which the 
vv -^ vv scattering can be coherent. This corresponds to the configuration in 
which the projectile and the scatterer switch places after scattering (Figs. (8b)). 
In contrast to the usual t-channel scattering, now the coherent condition allows 
the virtual Z to have an arbitrary momentum q^ ^ [p — p'"^) 

Analogous to the Bhabha scattering in QED, there is an s-channel graph for 
the ui? -^ vv scattering (Fig. 8d), in which the virtual Z carries a momentum 
(j9 + p')^ = —q^- Again, the momenta of the virtual Zs in Fig. 8b and Fig. 8d 
have the opposite sign, and thus their sum survives the cancellation that occurs in 
lowest order. 
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4.1. Ue + UeiUe) ^ Ue + UeiUe) 

For definiteness, we now consider the z/g^'e -^ Ve^e scattering. Of cosmological 
interest, we consider situations where |g| is of the order of a few MeV . The tree 
level matrix element of Fig. 8b is simply 

Mo = - (1) ' ^2 - 2 [^'-e iP)^o.Lu,^ (P)] [^^e (P)7"^^.. {P% (4.1) 

where p' and p refer to the momenta of the scatterer and the projectile, respectively. 
The over all minus-sign is due to the exchange of identical fermionic particles. 
Following the discussion of section 2, one finds that the index of refraction arising 
from neutrino-neutrino and neutrino-anti-neutrino scattering in a medium with 
equal amount of neutrinos and anti-neutrinos is 

/ i\ cos2 6'v^^ _,. ,. ^. 

where (n - l)i,^+e{e)^v,+e{e) is given by Eq. (2.10). 

Because of the absence of QED corrections, the calculation of one-loop radiative 
corrections to neutrino-neutrino scattering is considerably simpler. However, in the 
case of z/gC -^ ^e& scattering we are interested in radiative corrections to the charged 
current, whereas here the interest is shifted to the neutral current interaction. 
The result will therefore have a moderate dependence on the top quark mass. 
For convenience, we again orgainize our results into the self-energy, the vertex 
correction, and the box-diagram parts 

^^Weak ^ ^^^^^^ ^ dMvertex + ^Mbox- (4.3) 

The result for the self-energy and counterterm (Fig. 9) can be obtained from 
the analysis of neutral current radiative corrections of Marciano and Sirlin (Ref. 
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32) 

SMself = Mo 



Azziq^) - ReAzziMl) 25e j c^ \^(Azz{Ml) Aww{M^] 



fUr-lKRe 



q^ -Ml e ■ Vs2 7 V M| Mf 



w 
(4.4) 



where Mq is given by Eq. (4.1). 

Neglecting the neutrino and the electron mass, we find that the leading log 
vertex correction (Fig. 10) is 



5Mvertex = Mq—- 

4:71 



^M7-'«-'"i:i)-4(^)'"i^]' <^^^' 



Here, we have made the approximation indicated by Eq. (3.23) which is valid for 
our purpose. A more complete result without this constraint can be found in the 
appendix. In any case, one can see from Eq. (4.5) that the vertex correction to 
the neutrino index of refraction is completely negligible because 45^ — 1 is close to 
zero. 

Finally, there are a total of six box diagrams (Fig. 11). The sum of their 
leading log contributions is 

6MBo. = Mo-^(2-^+^-^^\n^). (4.6) 



47rs2 V 2c2 6 M^ q^ 

In reaching this result we have again made use of the approximation of Eq. (3.23). 
Also, non-coherent terms similar to those encountered in the above section have 
been excluded. It then follows that in the renormalization scheme discussed above 

1 + Ai? (g ) + ^^^2,2 6 1^^, 

(4.7) 



xU+ ' 



2 



M| 
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where AR'{q ) is the analogue of AR{q ) in neutrino-electron scattering 

Again only the real part will be kept. The over all normalization constant p'^'^) can 
be conveniently related to the constant Pj^'(j introduced in the study of ui -^ ui 



>;j^) = >;^) _ ^{^_\ [1 in„2 , io„4 



scattering (Refs. 32, 25) 

where 

tf} = 1 + - f 4 In c2 - 4. + S ■ -^'^ '^2^ , --. 






W 

8 2 
and i = my Ml. 

Again, the order a correction (the first term in Eq. (4.7)) cancels for the sum 
of neutrino-neutrino and neutrino-anti-neutrino scattering. Terms surviving the 
cancellation are grouped in the squared brackets. Among them, the first term 
is due to the expansion of the Z propagator, the second is due to the Z self- 
energy and the last term is the sum of vertex and box diagram corrections. The 
difference between p^^'^^' and p)^q arises because the three exchange box-diagrams 
(Fig. lid to Fig. llf) do not appear in a z/£ ^ vd, scattering. Numerically, this 
difference is very small, and can be ignored in our leading log approximation. The 
aforementioned top quark mass dependence is included in the parameter p\^Q (Eq. 
(4.10)). This gives an additional correction of the order +1% for a heavy top with 
mass 150 GeV <, mt <, 200 GeV . 
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In contrast to the VeG -^ ^e& scattering, no significant cancellation between the 
vertex correction and the box diagram correction occurs. As we already pointed 
out, this is because the vertex correction is suppressed by 4s^ — 1. As a result, 
within the region 1 MeV^ < |g2| < 400 MeV'^ , 6Mvertex + ^Mbox is essentially 
given by 5Mbox- Numerically, we find that it is about +13% to +10%. 

For the self-energy contributions, we again decompose them into the bosonic, 
hadronic, and leptonic parts 

Ai?'(g2) = Ai?'(^)(g2) + AR'^^\q^) + Ai?'W(g2). (4.II) 

The bosonic part does not have any large log terms, and hence its contribution is 
negligible. Explicitly, 



Ai?'(^)(g2) 



1 

- -, 9 , . /I 
a 



47rs2 



1/17C^ l + S^ 2 4 2\ 1 /^ , X(X-I) 

6 V 2 2c2 ) c^ J x{l-i)+i 



/17 S^ s2 r 2\ f , , C^+X{X-1) 

+ \ 7 - ^ - 2s^ + 5c2 / rfx In ■ 



V 2 2c4 c2 J J c2 



+ ^r— + -^ - 1 / dxx{x - 1) In ■ ^ ^ 



2 2c2 / / c 



L.2 



1 

2 



1 f.^f^2 , n .^. o^,.^' + e(l-^) 



H 2 I dx(x^ + {I - x)i - 2\ In 



X + ^(1 — x) 


(4.12) 
Numerically, we find that the terms in the square brackets add up to a value of the 
order of —1 and they are not sensitive to the choice of ^. Thus, AR'^ >(q'^) ^ 0.3%. 

The hadronic part suffers the same ambiguity as AR^ '{q^) induced by the 
complication of strong interaction. Within the constituent quark approximation, 

28 



we find 

3a >r^/l . o„ „ ..„o\„ f m^r-Mlx(l-x) 



Ai?'W(g2) = -^J2{l-'2s''CsfQf+4s^Q})Re f dxx{l-x)\n 



7 '''Z-' 



m 



2 



/ n -I 



Q'^X^l — x) I 



(4.13) 
where Qf and C^f are the charge and the third component of weak isospin of the 
quark / with a constituent mass rrif. Again, it is easy to show that the contribution 
due to a heavy quark is completely negligible in AR'^ ' {q'^) . Numerically, one finds 
from Eq. (4.13) that for \q^\ < m^, AR'^^\q^) is about +5%. 

The leptonic contribution, AR'^> can be obtained from Eq. (4.13) by an 
appropriate change of charge and weak isospin quantum numbers for the charged 
leptons and neutrinos, and then dividing the result by 3. The sum is over all the 
three families of leptons and neutrinos. In the region 1 MeV'^ < |g^| < 400 MeV^, 
we find that AR'^^\q^) is between +7% and +5%. 

Thus, the sum of all weak corrections plus a +1% correction if the top quark 
mass is within the region of 150 GeV to 200 GeV is quite sizable (~ +20%). The 
details are summarized in Table 1, in which one finds that for 1 MeV'^ < |g^| < 
400Mel^2 ^hg ^o^al correction is between +25% to +20%. 

The result for Uei^e -^ ^'e^'e can be obtained from that of i'e{p)ve{p') -^ ^e{p)ve{p')\ 
by a change of u{p')'jaLu{p') -^ —u{p')'yaRu{p') and replacing q^ by — g^. While 
the order a constant term will cancel for the sum of these two amplitudes, terms 
proportional to q^ /M"^ add up. Thus, the index of refraction due to Vg^eiyVf,) -^ 
z/eZ/e(z/e) Scattering in a medium which contains an equal amount of z/g and z/g is 

{n - 1).,.,(.,)^.,.,(.,) = -^G^p(^'^)iV,,po(Po)(l + 5...J. (4.14) 

In the leading log approximation, 5i,^i,^ is about +20% to +25% in the region 
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1 MeV < \q \ < 400 MeV . The details are summarized in Table 2. 

In the same region of q^ discussed above, we now have a slightly different 
hierarchy 

Mw,z ^ ml,r ^ \q\ (4.15a) 

However, the log factor \\i{M'y^/mf) {i = /i, r) is still rather large, and thus a 
leading log approximation calculation is still meaningful. In this case the function 
Ji defined in Eq. (3.1) can be approximately simplified to 

1 771 

/i(mj.„m^,Mw^,g2)~_ In— f. (4.156) 

6 M^ 

As a result, some of the log factors in 6Mvertex and 5Mbox in Eqs- (4-5) and (4.6) 
should be changed from \\i{q /M^r) to \\i{myMy^). The rest of the calculation 
remains the same. Places where such a change should be made are those in which 
diagrams contain charged leptons. 

For the vertex correction we find 



dMvertex = ^0^— 



-[-fV(l)-ln^U^r^£^ln^+^ln^^ 



(4.16) 

where 

ig'^ 1 
Mo = ■^j^^^—^[u^,{p')-faLu^Ap)]KAph'^Lu^Ap')]- (4-17) 

Eq. (4.16) reduces to Eq. (4.5) if one substitues mj by g^ and ignores the difference 
(valid in the leading log approximation) between M^j/ and Mz- 
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Accordingly, the box diagram contribution becomes 

It then follows from Eqs. (4.16) and (4.18) that for ug + ug ^ ug + ug {£ = fi, r) 
our leading log result is 

Mo + 5MW^^«^ = - ^v^G'^p(-'-)[M,,(p)7,Ltx,,(p)][M,,(p07"i^«..(p')] 



xn+^' 



M| 



Atts'^c'^ V 6 g2 g ^2 

(4.19) 

It bares a strong resemblance to Eq. (4.7) for z/g + z/g — * z/g + z/g scattering. 

Given the smallness of s^, the difference between Eqs. (4.7) and (4.19) is almost 
purely academic. Thus, qualitatively, we expect that the total correction for the 
three types of neutrino- neutrino scattering are approximately the same (~ +20%). 
The details are summarized in Table 2. 



5. Scattering Between Different Families 

In the absence of family mixing, lowest order interactions between different 
families such as z/^^t + e(e) -^ u^^r + e(e) and z/j + I'jipj) ^> I'i + J-'j{yj){i 7^ j) 
can only occur through the t-channel. Then, for the reasons explained above, the 
coherent amplitudes of these interactions will cancel completely, and they will not 
contribute to the neutrino index of refraction. The situation will change if one 
goes beyond the lowest order. In this case, radiative corrections become the only 
contribution. 
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5.1. iy^,T + e{e) ^ u^^r + e{e) 

It should be clear by now that in the absence of family mixing only those 
diagrams which are similar to the last four in Fig. 7 contribute to the neutrino 
index of refraction. The only difference is to change the z/e line by a u^^r line, and, 
accordingly, the internal electron line by a ;U or r line. The result is 



M.,,.+e^.,,.+e = ^O^^l^ + ^ 



■' 1 In — 75 — h-ln — ^-^ 



O , III' 

-\ — In — 
W 



Ml 3 g2 



(5.1) 
where the ellipses represent incoherent terms which are of no interest to us, and 

. 2 

Mo = -^[u,ip'h^Lueip%u,,M^''Lu,,M]. (5.2) 

In Eq. (5.1) terms proportional to q /Myy contribute to the neutrino index of 
refractions. By employing the by-now familiar method, we find 

(n - l)^^+eie)-^ue+e{e) = '^^^ ^NePQ{pQ)5j,,e, i = ^i,T. (5.3) 

In the region 1 MeV'^ < \q'^\ < 400 MeV"^, 6,y^e varies from about +11% to +10%, 
and 6,yT is between +7% to +6% (see Table 2). Notice also, A^e = 2A^jy (see Eq. 
(1.1)). Such a sizable correction in the diagonal scattering {ueC — > z/gc) is cancelled 
by the vertex correction. There is no vertex correction for scatterings between 
different families if the corresponding mixing is zero. 
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5.2. ui + yj^yj) —> ^1 + i'j{i'j){i ^ j) 

Leading log contributions are generated from box diagrams similar to the first 
three graphs in Fig. 11. A straightforward calculation shows that in the zero 
mixing limit 



M 






[Uu, {phaLUj,^ (p)] [Uy^ {p')-i'^LUy^ (p )] 



a 



X 



A'Ks'^C^ 



4^M2 Vl2 



— In- 



M 



+ -In- 



M 



m„ 



q' 



+ ..., 



(5.4) 



where the ellipses again refer to the incoherent terms and mi. is the heaviest 
charged lepton mass in question. It then follows that 



(n-r 



V^ + Vj(Vj)^Vi+l'j{l'-j) 



8^2 



Gf^Nj,^Po{Po)6j,^i 



(5.5) 



Numerically, we find that in the region 1 MeV^ < Ig^l < 400 MeV^, 5i^^i^^ is about 
+6% to +5% for the scattering z/g,^ + Ur ^ t'e,^ + t^r- The result increases about 
+ 1% for t'e + ^'/i ^ ^'e + '^fi- The details are summarized in Table 2. 



6. Neutrino-photon Scattering 

Since a photon is its own anti-particle, the cancellation discussed above does 
not apply to the neutrino-photon scattering. However, we will show below that 
the neutrino-photon forward scattering amplitude is zero at one-loop level in the 
standard model. 
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It is known that to the lowest nonvanishing order the diagonal effective 
neutrino-two photon interaction is given by 



L. 



eff 



auuF^^^F^f, + 



p + ia z/75 



vF^'^F, 



a/3; 



(6.1) 



where z/ is a neutrino field, F"'^ is the electromagnetic tensor with its dual F"^ = 
■^e^^" F(jX. In terms of the standard electroweak interactions, a and a' are zero at 
one-loop level if neutrinos are massless. Assuming neutrinos are massive, one then 
finds (Ref. 35) 

; Gpa 



■mj,f, a = 0, 



(6.2) 



where rrii, is a neutrino mass and / is a scalar function obtained from loop integrals 
with its real part given by 
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sm 



-1 



ki.k2/2mj 



AkiM 4:{ki.k2y 



2 2 ^^'■ 



1 ^^2 l+^/l-2mykl.k2 
l-^l-2mf/fci.fc2' 



if ki.k2 < 2mj 
if ki.k2 > 2mi 



(6.3) 



Here, ki and k2 are, respectively, the initial- and the final-photon momentum. 

The function / has a notable feature that it is zero for a forward scattering 
in which ki.k2 = 0. Thus, in the forward direction L^ff is zero at one-loop level. 
Notice also, the matrix element u{p)i^'y5u{p)iy is zero. We then conclude that the 
neutrino-photon scattering does not contribute to the real part of the neutrino 
index of refraction at one-loop level. Thus, the photon plasma is essentially trans- 
parent to neutrinos. 
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7. The Imaginary Part of the Index of Refraction 

While the real part of the index of refraction describes the coherent interference 
of propagating neutrinos, the imaginary part characterizes the incoherent depletion 
of neutrinos from their original coherent state. In normal situations, Re{n) ^ 
Im{n)^ because neutrinos only participate in weak interactions and hence the phase 
shift in the forward scattering amplitude is very small. In the MSW formalism, for 
example, Re{n) is of the order of Gi?, whereas Im{n) is of the order of Gpq . In a 
CP symmetric plasma such as the early universe, however, the order Gp of Re{n) 
cancels. As we have already learned, the leading terms of Re{n) are now only of 
the order of Gpq^ /Myi^. Since q^ <^ My^^ in this situation the imaginary and the 
real part of index of refraction become interestingly comparable. 

The simplest way of calculating the imaginary part of index of refraction is to 
employ the optical theorem 

Mn) = ^5^iV,a,-, (7.1) 



where Nj is the density of the jth scatterer and aj is the corresponding total cross 
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section. Table 3 lists all the relevant scattering cross sections normalized by 



.0 ^ ^a<^. (7^2) 

ovr 



The result given by neutrino-photon scattering is completely negligible. 
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8. Summary 

In this paper we have systematically evaluated the standard electroweak inter- 
action contributions to the neutrino index of refraction in the early universe. Of 
cosmological interest, we have concentrated on the period when the temperature 
of the universe was of the order of a few MeV, and the scatterers were photons, 
electrons, positrons and three types of neutrinos and anti-neutrinos. 

Assuming CP invariance, the number of particles and anti-particles is equal, 
and as a result the lowest order coherent forward scattering amplitudes completely 
cancel. Hence, the leading nonvanishing result obtained by the expansion of the 
gauge boson propagator is very small. In the absence of leptonic family mixing, 
which happens if Ue, v^ and Vr are degenerate in mass, even such a tiny result van- 
ishes for scatterings between different families. Although the cancellation does not 
occur for neutrino-photon scattering, we find that its coherent forward scattering 
amplitude is zero at one-loop level. 

The smallness of tree level result motivates us to investigate its radiative cor- 
rections. By employing the on-shell renormalization scheme, we have found that 
depending on the specific scattering process the leading log corrections are typi- 
cally of the order of 20% to 10% with the same sign. Radiative corrections become 
the dominant contribution for scatterings that involve different families if mixings 
in the lepton sector are very small. Our results for neutrino index of refraction are 
summarized in Table 2 and Table 4. The justification of our leading log approxi- 
mation relies on the hierarchy M^^ ^ ^ \q'^\^mf. 

Numerically, these corrections already become significant. They are about two 
orders of magnitude bigger than the original expectation (Ref. 11) of the order 
of a/TT. As summarized in Table 2 and Table 4, the radiative correction to the 
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neutrino index of refraction is about +20% and +50% for z/g and t'^,r, respectively. 
This generates a sizable effect in locating the exact boundaries from the nucleosyn- 
thesis constraint for the aforementioned Ve — ^s oscillation. A detailed discussion 
of the application of our results to cases of cosmological and astrophysical interest 
will be presented in a forthcoming paper. 

Besides the numerical significance, we have proposed a theoretical method 
to evaluate the neutrino index of refraction at finite temperature. This method 
makes the lowest order calculation much simpler than that derived from the finite 
temperature field theory, and more importantly allows us to evaluate higher order 
corrections in a straightforward way. The underlying physics also appears to be 
more intuitive and transparent. 

We also studied the question of the cancellation of infrared divergence in a 
coherent process. We have shown that such cancellations indeed take place in all 
order in renormalizable perturbation theory. As a consequence, the cancellation of 
infrared divergences will not spoil the coherent condition. 
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APPENDIX 

In this appendix we give some details of our calculation for the three- and 
four-point functions. As explained in the text, we have ignored all t-channel one- 
particle-reducible graphs. In the remaining graphs we keep all terms of order a 
and only leading log terms of order Q;(g^/M^). The order a terms, although 
they do not contribute to the neutrino index of refraction, provide a check on the 
calculation. Our calculation is carried out in the Feynman gauge. Due to the 
cancellation of particle-particle and particle-anti-particle scattering, graphs with 
scalar lines can be omitted in the leading log approximation. Consequently, only 
the leading log terms of our result are gauge invariant. 

Three-point Functions 

The three-point function correction to z/gC -^ Vet given by (JMge+e/ in Eq- 
(3.17) of the text is obtained by evaluating the one-loop vertex diagrams shown 
in Fig. 12. For convenience, we parameterize it as a sum of individual diagram 
contributions 

^5' r- / /^ r / M « 
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Similarly, for the vv -^ vv scattering (Fig. 13), we have 
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Box Diagrams 
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First, the photon box (Fig. 4c) result is 
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where A is a ficticious photon mass. 

In calculating the remaining weak interaction box diagrams, we frequently 
encounter a scalar integral 
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An explicit calculation shows 
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where the ellipses represent the terms which do not have a large log, and 
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They are related to the functions /i,2 (Eq. (3.1)) by 



Ko{m,m£,M,q^) = Ii{m,m£,M,q^) + l2{m,m£, M,q'^), 
M2 
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From the analyticity of scalar integrals, one can show that for our case only 
those graphs which have an s and/or -u-channel light intermediate state have a 
large log. For the UgC -^ v^e scattering (Fig. 7), we find 
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For the scattering z/z/ -^ vv (Fig. 11), we find 
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FIGURE CAPTIONS 

Fig.l. Tree level graphs for z/ge — > z/gC forward scattering. The coherent con- 
dition only allows the W in the w-channel (la) to carry a non-zero momentum. 

Fig. 2. Tree level graphs for z/gC -^ UeC forward scattering. Again, only the W 
in the s-channel (2a) is allowed to carry a non-zero momentum. 

Fig. 3. Feynman graphs for an electron forward scattering with an infrared 
finite but otherwise arbitrary potential (the blob) and its QED corrections. The 
ellipses represent other graphs which are infrared finite. 

Fig. 4. One-loop QED corrections to the -u-channel neutrino-electron forward 
scattering. The photon propagator in (4a) and (4b) is the massless photon prop- 
agator (the second term in Eq. (3.10)). The box diagram (4c) is evaluated with 
the full photon propagator. 

Fig. 5. One- loop self-energy and counterterm contributions to the w-channel 
neutrino-electron forward scattering. 

Fig. 6. Vertex and wave function renormalization corrections to the -u-channel 
neutrino-electron scattering. The photon propagators in (6a) and (6b) are the 
massive photon propagators (the first term in Eq. (3.10)). The detail of the blob 
is shown in Fig. 12. 

Fig. 7. Box diagram contributions to u^e -^ z/gC forward scattering. 

Fig. 8. Tree level diagrams for uu -^ vv and vv -^ vv forward scattering. 

Fig. 9. One-loop self-energy and counterterm contributions to vv -^ vv for- 
ward scattering. 

Fig. 10. Vertex corrections to the exchange channel in vv — > vv scattering. 
The detail of the blob is shown in Fig. 13. The ellipses represent the wave function 
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renormalization diagrams. 

Fig. 11. Box diagram contributions to vv -^ vv forward scattering. 

Fig. 12. One-loop vertex diagrams for Ve,e -^ v^^e forward scattering. Here the 
photon contribution is evaluated with the full propagator. 

Fig. 13. One-loop vertex diagrams for vv -^ vv forward scattering. 
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TABLE CAPTIONS 

Table 1. Radiative corrections to the neutrino index of refraction, assuming 
the medium is CP symmetric. Here |gp varies from 1 MeV'^ to 400 MeV"^. The 
definition of the various corrections is given in the text. A 1% correstion should be 
added to the diagonal uu -^ vv scattering if the top quark mass is within 150 GeV 
to 200 GeV . Results for scattering between different families in the small mixing 
limit are given in Eqs. (5.3) and (5.5) of the text, and are summarized in Table 2. 

Table 2. Summary of leading log results of radiative corrections. The interac- 
tion process also includes the corresponding scattering with the anti-scatterers. 

Table 3. The total scattering cross section for neutrinos with the various scat- 
terers in the early universe. The result is calculated in the center of mass systerm 
and normalized by o^ (Eq. (7.2)). The plus-sign refers to i = e and the minus-sign 
corresponds to i = fi, r. Also, i ^ j- 

Table 4. Neutrino index of refraction in the early universe. Here, no is given 
by no = —{2\/2/My^)G^N^{p'Q). The densities of the scatterers are normalized 
in terms of the photon density according to Eq. (1.1). The various radiative 
correction results, (5^, are summarized in Table 2. 
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Sourses 
hadronic self-energy:5%:5%:5% 
aic self-energy:4% to 3%:7% to 5%: 7% to 5% 
onic self-energy:negligible:negligible:negligible 
3rtex+box:< 2%:13% to 10%:12% to 10% 

QED:-l%: : 
Gal: 10% to 7%:25% to 20%: 24% to 20% 



z/ee(e) -^ Ueeie 



^eM'^e -^ '^eM'^e 



^IJ.,Tl^fl,T\l^fl,T ) 



^ ll,T^ ^,T\y ll,T } 



Table 1 
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correction: interaction: result ( 


q =\ MeV): result (|g| 


= 20 MeV) 


Su.u.-l^el^e — 


> z/ez/e:25%:20% 




^U^U^-l^^l^^Ji -^ 


z/^z/^: 24%:20% 




5u^j,^:UrUr - 


> Uriyr:2A%:20% 




5i,,v^,-VeV^ - 


-^ Ueiy^: 7%:6% 




^v^v^'-^e^T - 


-> z/e^r:6%:5% 




^v^Vr'-VpiVr - 


-^ z/^z/^:6%:5% 




^v.e-Ve(^ - 


>z/ee:10%:7% 




^v,,e-Vp.e. -^ 


z/^e:ll%:10% 




^Vre-Vre- - 


-^ z/^:7%:6% 





Table 2 
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interaction-channel: cross-section 


ViVi <-^ ee:8s^ ± As^ 


+ 1 


Vii>i '^ v^^y.\ 




Vie ^^ Vie:%s^ ± 6s^ 


+ i 


Vie ^^ Uie-.Ss^ ± 2s^ 


+ h 


UiUi ^ UiUf.e 




UiUj ^-^ UiUj-.S 




ViVi ^ UiUfA 




UiUj ^^ UiUj-.l 




z/(z/)7 ^^ z/(z/)7: negligible 



Table 3 
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real part of index of refraction: results 



(n - l).,/no:2(l + 5, J + \[c^ p^-'^-) {I + 5,^,J + 5,^,^ + 5,^,J 
(n - l)^,/no: 5,^e + \[c^ p^'''''\l + b^,v,) + 5.,^^ + ^v,Vr\ 



Table 4 
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